The contact conductance of a one-dimensional wire partly embedded in a 

superconductor 
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The conductance through a semi-infinite one-dimensional wire, partly embedded in a supercon- 
ducting bulk electrode, is studied. When the electron-electron interactions within the wire are 
strongly repulsive, the wire efltectively decouples from the superconductor. If they are moderately 
or weakly repulsive, the proximity of the superconductor induces superconducting order in the seg- 
ment of the wire embedded in it. In this case it is shown that the conductance exhibits a crossover 
from conductive to insulating behavior as the temperature is lowered down. The characteristic 
crossover temperature of this transition has a stretched exponential dependence on the length of the 
part of the wire embedded in the superconductor. The amount of this stretch is determined by the 
nature of the electron interactions within the wire. 



I. INTRODUCTION 



In recent years one dimensional interacting electron 
systems have attracted a large amount of attention. Part 
of the interest in these systems lies in the fact that 
electron-electron interactions in 1-D wires, even when 
weak, cannot be considered perturbatively. A ques- 
tion of practical importance, dealt with by a number 
of authors^^, is that of the contact conductance of a 
one dimensional system connected to an external elec- 
trode. Most works on the subject picture the junction 
as a one dimensional wire connected to the electrode at 
a point. Here it has been found that electron-electron 
interactions strongly influence the conductance through 
the junction. In particular, repulsive interactions in the 
wire drive the system to be insulating at low enough tem- 
peratures, unless the contact is perfectly clean. Attrac- 
tive interactions, on the other hand, mask obstructions 
at the interface between the two systems. These results 
do not qualitatively change whether the electrode is su- 
perconducting or metallic. 

In this paper we explore the behavior of a junction 
between a 1-D wire and a superconducting electrode of 
different geometry, specifically, the situation where the 
wire is embedded some distance into the electrode, as 
illustrated in Fig. [TJ In this case one expects that super- 
conducting order is induced in the part of the wire which 
is embedded in the superconductor, enhancing the con- 
ductance of the junction. On the other hand, for point 
contacts with even a small amount of normal reflection 
(and repulsive interactions in the wire), it has been shown 
that superconducting order suppresses the conductivity 
of the junction^'^. The main goal of this work is to clarify 
how these two competing effects determine the behavior 
of the junction as a function of the temperature (or the 
applied voltage), and length of the embedded wire. 

As expected, provided interactions within the wire are 
not too strong, superconducting order is induced in the 
part of the wire embedded in the bulk, and an effective 
gap is formed in the wire whose value is determined by 
the tunneling rate between the superconductor and the 
wire, as well as the nature of the electron interactions. 
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Figure 1: A schematic representation of a 1-D wire in con- 
tact with a superconductor: (a) A point contact, (b) A wire 
embedded a length L into the bulk 



Yet, in spite of this proximity effect, the junction conduc- 
tance exhibits a crossover from a conductive state, gov- 
erned by Andreev scattering, to an insulating behavior at 
low enough temperatures, characterized by a power law 
dependence on the temperature. This behavior results 
from the finite amount of normal backscattering within 
the wire. We found that the characteristic crossover tem- 
perature between the two regimes has a stretched expo- 
nential dependence on the length of the junction, and the 
amount of stretch is determined by the strength of the 
electron interactions in the wire. 

The article is organized as follows: In the next section 
the model whereby a 1-D wire is embedded infinitely deep 
into a bulk BCS superconductor is introduced. By inte- 
grating out the superconductor degrees of freedom the 
effective action of the embedded wire is obtained. In 
Sec. II and the Appendix, renormalization group (RG) 
techniques are employed in order to characterize the be- 
havior of a 1-D wire embedded in a BCS superconduc- 
tor. This RG flow allows one to deduce the low energy 
properties of the finite part of the wire embedded in the 
superconductor. In Sec. Ill the conductance of the wire- 
superconductor junction is evaluated. Finally, the results 
are summarized and discussed in Sec. IV. 



II. MODEL 

We consider a single 1-D wire with interacting electrons 
(including backscattering), embedded inside a standard 
BSC superconductor (SC). The action of the system is a 
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sum of three contributions: 

S = Ssc {(p, ^) + Sw ('0, -0) + St (4), ip; (p, if) (1) 

where Ssc is the action of the superconductor, Sw is the 
action of the wire, and St describes the tunnehng between 
the two systems. is the electronic field operator in the 
wire, (fia- is the field operator in the SC, and a denotes 
the spin index. 

The action of the SC, modeled by the standard BCS 
hamiltonian with a constant pairing amplitude A^c, has 
the form 



dr + Hq Asc 
A* dr - Ho 



(2) 



where the vector ^ = (x, y, z, r) contains three space co- 
ordinates and an imaginary time, <i>^ = (f^) is the SC 

electron field in Nambu notation, and Hq = ^^^V^ — /i 
is the free Fermi gas Hamiltonian, with m and fi as the 
electron mass and chemical potential, respectively. 

The 1-D wire is modeled as a Bosonized Luttinger liq- 
uid. Since the literature on this subject is quite extensive 



(ref. [13] and therein), we will only point out the more 
relevant details to the topic at hand. 

We represent the electronic fieldai^, 



V'r,o-(a;) 



V2; 



1 Jrkjx^-iy/^[r(l>a{x)~ea{x)+ij{r<l>4x)~e4x))] 

(3) 

in terms of boson charge and spin density fields: dxtj^^ix), 
and their conjugates, 6^{x), where v = c/s denotes the 
charge and spin sectors, respectively. In the above for- 
mula, kf \s the Fermi wave number, a is the short dis- 
tance cutoff, r is the chiral index representing the right 
{r ~ 1) and left (r = —1) moving part of the electronic 
field, and a — 1 for ] spin while cr = — 1 for | spin. 

In terms of the boson fields, the action of the wire 
becomes a sum of three contributions: 



Sw = Sc{Oc) + Ss{4is) + Sbs{4>s), 



(4) 



where 



Sa{ec) 



dx 



drUrKr 



(5) 



describes the charge sector (after the dependence on the 
(j)c field has been integrated oui*i^). A similar expression 
describes the spin sector, Ss{4>s)^ with the subscript c re- 
placed by s and Kc replaced by Kj^. Here Kc^ Kg, Uc, 
and Us are model specific parameters describing the in- 
teraction strength {K) and the mode velocity (u) of the 
charge and the spin fields. For the noninteracting case 
Kc = Ks = 1 and Uc = Us = Vf, where Vf is the Fermi 
velocity of the wire. For an interacting system, values 
of Kc > 1 and Kg < 1 correspond to attractive inter- 
actions, while Kc < 1, Ks > 1 correspond to repulsive 
ones. Generally in an interacting system the velocity of 
the two modes differ, Uc ^ Us- 

The third term of the wire's action, Sbs, describes 
backscattering of two electrons with opposite spins. 
Namely a collision which effectively results in a spin flip 
between the right and left moving parts of the electronic 
field (~ a'^R,'y''l'R -u'4^L,-iy)- This term has the form 



Sbs — 



2.9 



(27ra)' 



dxdr cos 



(V8^0,(x,r)) , (6) 



where g is the backscattering coupling constant. 

Finally, the tvmneling between the superconductor and 



the wire is described by 

St = j d^S,' j dxdr J2 ^a{e)t{e;x,T)i;„{x,T) 



H.C.. 



(7) 



where — V'r.o-, and t{S/ ; x, r) is the tunneling ma- 
trix element. In the simplest case this tunneling is instan- 
taneous, homogeneous in space, short distant and SU(2) 
spin symmetric. Under these assumptions it takes the 
form 



t{^'; x, t) = i{x)6{x' - x)S{y')S{z')d{T' - r) 
where t{x) has the characteristic function form: 



iix) 



to for < a; < L 
otherwise 



(8) 



(9) 



where L is the length of the part of the wire embedded 
in the superconductor (see fig. [T](b)). 

The effective action of a wire embedded in a super- 
conductor is obtained by tracing out the superconductor 
degrees of freedom. 



(10) 
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where S is given by ((T]). Since the action is quadratic in 
the superconductor field, this integration is straightfor- 
ward. The result may be written as a sum of three terms, 
Seff = Sw + Sph + Spp, where Sph, and Spp represent 
contributions associated with Giever (particle-hole) and 
Andreev (particle-particle) tunneling, respectively. At 
energies above the superconductor gap A^c the particle- 
hole term is dominant and it's contribution, on integrat- 
ing it out, will be to renormalize the chemical potential. 
On the other hand, in the low temperature regime (well 
below the superconductor gap Ag^) the particle-hole term 
vanishes and the main contribution comes from the An- 
dreev tunneling term, Spp. In this limit, after averaging 
over the rapid spatial oscillations, the tunneling becomes 
local in space and time. Expressing it in terms of the 
bosonic fields we thus have 



^eff - S 



W 



S 



pp 



(11) 



where Sw is given by (jH), and 

Spp = -— J dxdTCOS^y2TT(j}s{x,T)jcos (y2n0c{x,Tjj , 

(12) 

where A ~ (^toNojjj ■ Here A^'o is the normal density of 

states of the SC at the Fermi level, and p/ is the Fermi 
momentum. 

In terms of the original fermionic fields, the tunnel- 
ing term has the form of the regular pairing term in the 
standard BCS theory: 

Spp = aJ da;dr(-0|(x,T)V;^(a;,r) + Vl(a;,T)V'T(a;,r)) 

In the absence of electron-electron interactions, this 
term, along with the free quadratic kinetic term of the 
model can be diagonalized by the standard Bogolubov 
transformation. The excitation spectrum of this system 
will be gapped with an energy of A. 



III. THE RENORMALIZATION GROUP FLOW 
EQUATIONS 

The effective action dTTT) describes the physics of the 
junction for a temperature up to the order of the super- 
conductor gap, Asc (which will subsequently provide the 
high energy cut-off of our system). 

In order to describe the behavior of the system at much 
lower energy scales, and take into account the electron- 
electron repulsive interactions, we shall employ a real 
space RG approach, following Giamarchi & Schulz^-?;-— . 
As usual in these cases the RG procedure manifest it- 
self in a flow of the coupling constants of the problem 
as the ultraviolet cut off is reduced from 1/a to 1/a'. 
In our problem these coupling constants are: The in- 
teraction strengths, K^; The mode velocities, Wi,; The 
dimensionless backscattering constant, v = -2—; And the 
dimensionless pair tunneling strength A = A;^. The 
flow equations of these coupling constants (see Appendix 
for the details of the derivation) are: 



^ = xJ^A 
al V Uc 



dl 
dA 



-K^ f A,,A2 + ^ 



(13) 
(14) 



A{2--{Ks + K-'+y)) (15) 



dy 
dl 

dUr 



y(2 - 2K,) - 2A,A2 

2 



UeK-'W, ^A 
dl \uc 



dus 
~dl 



= UsKsWsA'' 



(16) 
(17) 
(18) 



where dl — d log a is the dimensionless change in the 
ultraviolet cutoff, 



^c(.) - ^ fj dv { cosHv) +{^) sm'{<p) ) (19) 



2ir / \ 2 



I^c(.) = d^cos(2(p) X |^cos2((^)+ (^^^ sm\^)j (20) 

and gc ~ while gs = Kg. 



We shall restrict our analysis to the spin symmetric 
case (i.e. the situation where the interactions between 
electrons with parallel and opposite spins are identical). 
In this casei^ Kg ~ 1 -|- 1 , and equations and for 
the dimensionless backscattering y and spin interaction 



strength Kg reduce to: 

^^V-2A.A^ (21) 
thus maintaining spin invariance. 
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A. Analysis of the Flow Equations 

The above RG equations imply that the charge and 
spin velocities, u^. and u^^ renormalize towards each 
other. This is a consequence of the correlations between 
spin and charge excitations generated by the proximity 
effect. From Eqs. P7|) . (ITOl) . and (^0)1 one can observe 
that the signs of Wc and Wg are determined by the ratio 
Uc/us in such a way that the velocities approach each 
other. 

Another consequence of the RG equations. ([T3|) and 
(fT4|) . is that Kc can only grow, while Ks can only be 
reduced. This flow^stem from the finite value of the tunj 
neling parameter A (and y). Thus the RG behavior of A 
controls the behavior of the system. 

The separatix between the regions where A is relevant 
or irrelevant can be obtained numerically for a given set 
of bare parameters. It is clear that for initially repulsive 
interactions (ifi"^ > 1, > 0, K^P'^ < 1), A will be 
relevant if: 

4-(ifi") + (i^(°))"'+y("))>0, (22) 

since the flow equations can only drive K^/s to be more 

" attractive" . In Fig. [2]we depict the RG flow in the A-Kc 
plane, in the spin symmetric case, for several initial val- 
ues of in the range between 0.32-0.325. Fig. [3] shows 
the separatix between relevant and irrelevant tunneling 
for initial bare (repulsive) values of Kc and y. One can 
observe that Kc must be smaller than ^ for the tunneling 
to be irrelevant. 
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Figure 2: Plots of the flow of A as a function of Kc for initial 
values of K^P^ between 0.32-0.325^^ The transition between 
irrelevant (A —> 0) and relevant (A oo) tunneling occurs 
at ~ 0.3223. The initial values of A'°) and y'-°'> are 0.05 
and 0.1 respectively 

In the case where A is irrelevant, the wire effec- 
tively decouples from the superconductor. The pair tun- 
neling between the two systems is suppressed and the 
superconductor-wire junction becomes insulating. Notice 
however that this behavior takes place at very strong re- 
pulsive interactions, Kc < ^, where the system tends to 
Wigner crystallize^. 
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Figure 3: The curve describes the separatix between relevant 
tunneling, where the system flows to a superconducting state, 
and irrelevant tunneling, where the system flows to a decou- 
pled wire, as a function of the bare values of Kc and y 



In the situation where A is relevant, the system flows 
to a singlet superconducting state. The interaction pa- 
rameters become attractive and spin-charge separation is 
no longer valid. Yet, before the interaction parameters 
obtain their asymptotic values (i.e Kc oo, Ks — *■ 
and y — )■ — oo which correspond to "infinitely" attractive 
interactions) , the RG equations [T^ will cease to be 
valid since the small parameters of our perturbation the- 
ory, A and y, will fiow to the strong coupling regime. 



B. Length Dependance of the Effective Gap 

In what follows we consider a weakly interacting Lut- 
tinger liquid, where the bare interaction parameters are 
close to unity. Moreover, in order to ensure that the per- 
turbative RG equations p^[T5|) remain valid we shall as- 
sume that L, the embedding length, is of order or smaller 
then the bare superconducting correlation length, Vf/A. 
Since L serves as an infrared cutoff for the RG fiow, this 
condition ensures that the this fiow is confined to the 
perturbative regime. We shall also assume that the tem- 
perature is much smaller than A, and approximate the 
velocities of the spin and charge sectors by their asymp- 
totic renormalized values: Uc = Ug = Vf. 

In the limit of weak interactions the RG equations for 
Kg/c a-re of second order in the perturbation parameters 
(A < l,y < 1), while the equation for A (Eq. [TS]) is 
of first order (neglecting the y dependance, which is of 
second order as well). The equation for the coupling y 
(Eq. [16]) is also second order in small parameters, at least 
if we consider the spin symmetric case . Thus one can 
assume Kc, Ks, and y to be approximately constants, 
and consider the simplified equation for A: 
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where 



7 ^ 



k: 



(24) 



Integrating the above equations from I = Q to I = 
log(-L/a), and using the relation A = Aa/z;^, we obtain 
the renormalized value of the gap of the wire embedded 
in the superconductor: 



A.,,.A(^ 



7-1 



(25) 



In particular, repulsive interactions [Kc < 1 and Kg > 1, 
and therefore 7 — 1 < 0) reduce the effective gap in the 
wire. 

This renormalization of the gap implies that the effec- 
tive correlation length. 



ieff 



1-7 



(26) 



is larger than the bare correlation length ^ = Vf/A. 

The above results are valid as long as the RG flow stays 
within the perturbative regime, namely A < 1. This con- 
dition implies that L should be shorter than ^{^/a)~~^ . 



IV. THE WIRE-SUPERCONDUCTOR 
JUNCTION 

At this stage of the RG procedure (pursued up to the 
scale L) the junction between the wire and the supercon- 
ductor may be assumed to be point-like. Thus if the tem- 
perature (or the applied voltage) is smaller than hvf/L 
one may continue to integrate out the high energy de- 
grees of freedom in the part of the wire which is not em- 
bedded into the superconductor down to the relevant en- 
ergy scale. This may be achieved following the procedure 
described in the literature^. The important ingredient, 
now, is the magnitude of normal back scattering from 



the junction. The latter is of order of 



roe 



where rg is system specific reflection amplitude in the 
absence of superconductivity, while S^^ff is the correla- 
tion length (|26p within the part of the wire embedded in 
the superconductor. This behavior of rjv results form the 
fact that only the charge that is not converted to the con- 
densate backscatters from the edge of the wire embedded 
in the superconductor. According tc^ the normal current 
reduces exponentially with the distance on a length scale 
of ^ , which gives the above estimate for r^r. Now, since 
^e// depends on L, the magnitude of the normal reflec- 
tion has a stretched exponential dependence: 



TN ^ roe 



(27) 



where a = 2a^^^ and the amount of stretch, 7, is dic- 
tated by the nature of the electron-electron interactions 
within the wire, as follows from 



The behavior of the backscattering, clearly, manifests 
itself in the conductance of the junction. In the limit of 
sufficiently high temperatures, one obtains^: 



where G%,g 



a 



9^0 



NS 



6G, 



(28) 



27Tn 



represents the conductance 



of an ideal junction where only perfect Andreev reflec- 
tions take place, and 



6G 



(29) 



The above formula holds for the range of temperature 
where SG <C G'j^g, since SG cannot be larger than G'j^g. 
Nevertheless, it has been shown^^, that any scatterer, 
at a point contact between a wire with repulsive inter- 
actions and a SC, will eventually drive the conductance 
to zero as the temperature is lowered. Since any finite 
length junction will have some backscattering, the con- 
ductance should drop to zero for low enough tempera- 
tures, as illustrated in Fig. 4. The crossover tempera- 
ture, T* from the conductive and the insulating behavior 
of the junction depends on the length of the wire embed- 
ded in the superconductor and its scaling behavior may 
be deduced from ([M]) and ([29ll : 



(30) 



This result implies that the temperature scale at which 
the effects of backscattering becomes substantial reduce 
as a stretched exponent with the length of the junc- 
tion, and the stretch is determined by the interactions, 
through the parameter 7. 



V. SUMMARY AND CONCLUSIONS 

In this work we studied a junction of a 1-D wire em- 
bedded a certain length, L, into a bulk superconductor. 
We first characterized the nature of the contact between 
the 1-D wire and the superconductor using a real space 
RG scheme. We found that repulsive interactions in the 
wire compete against the superconducting order being 
imposed by the bulk superconductor. The system can 
flow to either of two phases, depending on the nature of 
the interactions. When the interactions are strongly re- 
pulsive the tunneling between the two systems becomes 
irrelevant and the wire essentially decouples from the 
bulk superconductor. For moderate repulsive interac- 
tions, and for attractive interactions, tunneling is rel- 
evant, and the bulk superconductor induces supercon- 
ducting order in the wire. The gap opened in the wire 
depend on the tunneling strength, and electron-electron 
interactions modify its nominal value. 

The finite length of the part of the wire embedded in 
the superconductor, L, implies that the RG flow, in gen- 
eral, does not reach its asymptotic (non-perturbative) 
limit. Thus L introduces itself in the behavior of the ef- 
fective gap, and the effective correlation length, in the 
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Figure 4: A schematic graph of the conductance as a function 
of temperature for repulsive interactions in the wire. The de- 
cay of the conductance depends on the length of the junction 
as a stretch exponent. The longer the junction, the lower the 
temperature at which the decay sets in. In the limit of weak 
backscattering, represented by the solid curve, we can esti- 
mate the temperature dependance of the conductance through 
Eq. I|29|) . In the zero temperature limit, represented by the 
dashed curve, any initial finite backscattering will eventually 
drive the conductance to zero. The crossover temperature, 
Eq. (|30|l between the two behaviors is estimated to scale as a 
stretched exponent in the length of the junction. 

wire. This, in turn, dictates a stretched exponential be- 
havior of the normal reflection from the junction as func- 
tion of L. For energy scales (temperature or voltages) 
beneath the effective gap, we described the qualitative 
picture of the conductance as a function of temperature 
and length of the embedded segment of wire, see Fig. IH 

Several simplifications have been used for our analysis: 
One is that our model treats a semi-infinite wire with 
a single junction, while in practical situations a finite 
wire is usually connected to two reservoirs. This ideal- 
ization holds as long as the segment of the wire outside 
the superconductor is long enough compared to hvf/T, 
where Vf is the Fermi velocity and T is the temperature. 
Additional simplification is the assumption that the em- 
bedding of the wire into the bulk superconductor does 
not introduce inhomogeneities, i.e. the wire can be still 
considered to be clean, and that the tunneling to the 
superconductor is homogenous along the wire. This ap- 
proximation holds when the transport mean free path in 
the wire is longer than L. In the opposite limit one ex- 
pects a different behavior of the proximity effect which 
will change the conductance of the system. 
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VI. APPENDIX 

The mathematical formulation used in this work fol- 
lows closely the real space RG procedure used by Gia- 
marchi & Schulai^i^^: In this procedure one evaluates a 
correlation function in the wire, of the form (the time 
ordering symbol is suppressed): 

RUxa,Ta;xt,n) = ^env^M-..-»)-^(-..-.))^, (31) 

where tp can symbolize any of the boson fields and 7 is 
some constant. For the following discussion it will suf- 
fice to examine only one of the sectors, for instance the 
spin sector. The same considerations can be carried on 
straightforwardly for the charge sector. 

The fact that the relevant boson fields in the spin sector 
is 0s leads naturally to the evaluation of the correlation 
f mictions: 

Unfortunately this correlation function cannot be calcu- 
lated exactly using the complete effective action pTt . 
Though, if the tunneling and backscattering parameters, 
{t and g respectively), are small then it may be computed 
perturbatively. To second order in these parameters this 
function is found to be: 

i?s(f,,,) = ,-y'{K"^^irX.} + Dffs,n^i^,^^^,^))^ (33) 

where ra.b — ^ 'r'b, and (pp^s is the angle between the 
vector f = (x, u^t) and the x axis. The function is 
(at zero temperature)^'^: 

1 f x^ + iu,\r\+af \ 
Fs{x,T) = -\n\ — I. (34) 

Apart from the term proportional to sin^ {(pp^ ^^^g), this 
functional form is identical to the free correlation func- 
tion: 

i?(0)(r- - n) = e-^^^=^=(^^-^~^), (35) 

but with an effective Luttinger interaction constant K^^-f 
modified by the perturbations: 



a \a/ Z a \a 



r 



Here A = is the dimensionless tunneling parameter, y — is the dimensionless spin backscattering param- 
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eter, and Xs is a geometrical term given by: 



C0S^((y5)- 



sin^((y5) 



(37) 



Integrating out high energy degrees of freedom, near 
the ultraviolet cutoff, corresponds to integrating out a 
"small ring" between a a' = a + da, where a is the 
small distances parameter of the model. After integration 
and rescaling, an infinitesimal change is generated in the 
expressions for K^-f^ . In order to keep K^-l'-f constant 
with the reduction of the cutoff, it is required that the 
bare parameters change. For instance, we find that: 

K,{a')=Ks{ayK^,{a) ^\a)Xs{a) + ^^^ ^ (38) 
which generates an (exact) differential equation for Ks'. 

^ = -^,?(0(^.A2(0 + ^) (39) 

(here dl = 

In a similar fashion one obtains differential equations 
for the parameters A and y: 



dA^ 
dl 

dl 



y\4-AKs). 



(40) 
(41) 



The sin^((/?f7^ ^^^s) contribution to the correlation func- 
tion arises from the fact that the tunneling perturbation 
couples the spin and charge sectors, which were uncou- 
pled without this term. Mathematically, this term char- 
acterizes the anisotropy between the space (x) and time 
(mst) directions. It's pre-factor D'^^^ is given by: 



Dl^f ^Ds + KlK^wJ -(-)' , (42) 

where Ws is another geometric term factor: 



1 f^^ / ru \ ^ 



27r 



dipcos(2ip) cos (1^9) -I- 



sin'((p) 



Applying the same renormalization scheme to D'^^^ , 
will generate the flow of the parameter D^: 



dPs 
dl 



= KUl)WM)A^{l). 



(43) 



It should be noted that Dg is initially zero but is gener- 
ated under renormalization. 

The parameter Ds controls the renormalization of the 
velocity parameter. As long as the space and time di- 
rections are isotropic the velocity parameter does not 
flow under renormalization, but they should flow in the 
anisotropic case. Indeed, assuming that initially the cor- 
relation function is described by the function KgFg(fafi), 
then a small change of dug will generate the term: 
s\v?{ipp^ 6,s) ■ dug . Up to a factor, this is exactly the 
anisotropy term. Therefore, the renormalization of Dg 
is equivalent to that of the velocity Ug by the following 
relation (e.g. Eq. ES]):^ 

The flow equations for the charge sector can be ob- 
tained by an identical procedure. The equations obtained 



^ = = x,m ['^m) . (44) 

Finally these equations cannot be exactly correct, since 
they do not maintain the spin invariance SU(2) of a model 
that was spin invariant to begin with. (A spin symmetric 
model is one where the interactions between electrons 
of opposite and parallel spin are identical). Since the 
perturbations do not break this symmetry, something in 
the above result is insufficient. Indeed, it turns out that 
the remedy for this problem lies in the inclusion of the 
third order terms of perturbation theory. This correction 
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2X,A2. 
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